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ABSTRACT 

We present two-dimensional numerical simulations of the evolution of a low-mass star moving super- 
sonically through its surrounding interstellar medium (ISM). We show that the ejecta of a moving star 
with a systemic velocity of 20 kms -1 will interact with the ISM and will form bow-shock structures 
qualitatively similar to what is observed. We find that, due to ram-pressure stripping, most of the mass 
ejected during the AGB phase is left downstream of the moving star. As a consequence, the formation 
of the PN is highly influenced, even at the low relative velocity of the star. The models are based on 
the predictions of stellar evolution calculations. Therefore, the density and velocity of the AGB and 
post-AGB winds are time dependent and give rise to the formation of shock regions inside the cavity 
formed by the previous winds. As a result, the stand-off distance is also time dependent and cannot be 
determined by simple analytical arguments. 

Subject headings: hydrodynamics-ISM: structure-ISM: jets and outflows-planetary nebulae: 
general-stars: AGB and post-AGB-stars: winds, outflows 
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ties, large relative velocities or both were needed to explain 
PN asymmetries due to interaction with the ISM. 

This is the scenario commonly accepted nowadays and 
as a result, the search for interaction of PNs with the ISM 
has been mostly restricted to old PNs with large angu- 
lar extent (Tweedy & Kwitter 1996; Xilouris et al. 1996). 
Evidence for the interaction has been also reported in the 
study of some individual objects (e.g. Borkowski, Tsve- 
tanov, & Harrington 1993; Tweedy, Martos, & Noriega- 
Crespo 1995; Soker & Zucker 1997) and the observations 
have been interpreted in the light of the available models. 
As a consequence high spatial velocities were invoked to 
explain the morphologies of PN at large distances from the 
Galactic plane. 

The problem is that in the available models, the PN-ISM 
interaction was studied by considering the relative move- 
ment when nebular shell was already formed. In this letter, 
we demonstrate that when the PN formation is followed 
by considering the evolution of the stellar wind from the 
AGB, high velocities, high ISM densities or the presence of 
a magnetic field are not necessary to explain the observed 
asymmetries in the external shells of PNs. Moreover, in 
contrast to what is commonly believed, the interaction is 
observable even during the very early stages of the PN 
phase. 

2. THE NUMERICAL METHOD 

We have carried out the numerical simulations with the 
fluid solver ZEUS-3D (Stone & Norman 1992a,b; Stone, 
Mihalas, & Norman 1992), developed by M. L. Norman 
and the Laboratory for Computational Astrophysics. The 
computations have been performed using a 2D spherical 
polar grid with the angular coordinate ranging from to 
180° and a physical radial extension of 2 pc which gives a 
total grid size of 4x2 pc. We have adopted a resolution of 
400 x 360 zones in the radial and angular coordinates of the 
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1. INTRODUCTION 

The process or processes that may cause the departure 
of spherical symmetry in PNs have been the subject of 
many studies e.g. the hydrodynamic collimation of the 
fast wind by an equatorial density enhancement (Frank, & 
Mellema 1994), the rotation of a single star with a magne- 
tized wind (Rozyczka & Franco 1996; Garci'a-Segura 1997; 
Garcfa-Segura et al. 1999) and the evolution of a PN in 
a close binary system (Soker 1996). These mechanisms 
are able to account for the formation of elliptical, point- 
symmetric and bipolar PNs. However, several PNs have 
been observed in which mainly the external shell shows 
departure from sphericity (e.g. MA 3, PM 1-295, M 2- 
40). In these PNs, the presence of asymmetries cannot be 
explained by the physical processes related to the central 
star (CS), such as rotation, magnetic fields or binarity. In- 
stead, the interaction with the ISM was proposed to be the 
cause of the observed asymmetries. 

The interaction of a PN with the ISM was first sug- 
gested by Gurzadyan (1969) and the first theoretical study 
was made by Smith (1976); he assumed a thin shell ap- 
proximation and the 'snow-plough' momentum-conserving 
model of Oort (1951). Later on, Isaacmann (1979) used 
the same approximation, but considering higher velocities 
and ISM densities. Both theoretical works arrived at a 
similar conclusion, namely that the nebula fades away be- 
fore any disruption of the nebular shell becomes notice- 
able. Borkowski, Sarazin, & Soker (1990) analytically, and 
Soker, Borkowski, & Sarazin (1991) numerically confirmed 
the thin-shell approximation originally developed by Smith 
(1976), concluding that the interaction becomes dominant 
when the density of the nebular shell drops below a criti- 
cal limit, due to the nebular expansion. As a consequence, 
the interaction should be easily observable during the late 
stages of the nebular evolution. Moreover, high ISM densi- 
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grid respectively. We consider the photoionization of the 
gas during the PN phase by using a simple approximation 
to derive the location of the ionization front for arbitrary 
density distributions (see Bodenheimer, Tenorio-Tagle, & 
Yorkc 1979; Franco, Tenorio-Tagle, & Bodenheimer 1990). 
This is done by assuming that ionization equilibrium holds 
at all times and that the gas is fully ionized inside the 
H II region. The models include the Raymond & Smith 
(1977) cooling curve above 10 4 K. For temperatures be- 
low 10 K, the unperturbed gas is treated adiabatically 
but the shocked gas region is allowed to cool down with 
the radiative cooling curves given by Dalgarno & McCray 
(1972) and MacDonald & Bailey (1981). The photoionized 
gas is always kept at 10 4 K, so no cooling curve is applied 
to the H II region unless there is a shock. 

2.1. The Boundary Conditions 

We assume that the ISM moves relative to the star per- 
pendicular to the line of sight by fixing the position of the 
star at the center of the grid and allowing the ISM to flow 
into the grid at the outer boundary from to 90°. On 
the other hand, from 90 to 180° we set an outflow bound- 
ary condition. The temporal evolution of the stellar wind 
for a 1 M Q star during the AGB phase (taken from Vas- 
siliadis & Wood 1993) and the PN stage (Vassiliadis & 
Wood 1994) has been set within a small (five radial zones) 
spherical region centered on the symmetry axis, where re- 
flecting boundary conditions are used. The gas evolution 
in a static configuration as well as the stellar wind inputs 
for the models are described in detail in Villaver, Garcfa- 
Segura, & Manchado (2002a) and Villaver, Manchado & 
Garcfa-Segura (2002b). 

In order to study the interaction process we have 
adopted a purely hydrodynamical scheme, in which we 
assume that the ISM pressure is simply the standard gas 
thermal pressure, P = nKT (where n is the number of par- 
ticles per unit volume, K is the Boltzman's constant and 
T is the gas temperature). Since the strength of the inter- 
action depends linearly on the ISM density and quadrati- 
cally on the relative velocity of the star, we have adopted 
two ISM densities and a very conservative value for the 
star's velocity, 20 kms -1 . The ISM is assumed to be ho- 
mogeneous with the characteristics of the warm neutral 
medium (WNM) which is the main observed constituent 
of the ISM, with temperatures estimated to be in the range 
5000-8000 K (Kulkarni & Heiles 1988) and neutral atomic 
hydrogen densities between 0.1 and 1 cm -3 (Burton 1988). 
We use an ISM temperature of 6000 K and densities of 0.1 
and 0.01 cm -3 , the latter to study the interaction with a 
very low density ISM. 

3. RESULTS 

3.1. The Evolution during the AGB through a Low 
Density Medium 

The logarithm of the gas density during the AGB 
phase is shown in Figure 1 for a star moving with 20 
kms -1 through a low density ISM (n Q = 0.01 cm -3 ). Fig- 
ure 1 shows only half of the r — 9 plane, the star is fixed in 
the grid and the ISM flows in from the top to the bottom. 

During the very early stages of the evolution of an AGB 
star a free-streaming steady stellar wind with a velocity of 
2 kms -1 and mass-loss rate of M = 10 -8 M Q yr -1 reaches 



ram pressure equilibrium with the ISM at a stand-off dis- 
tance from the star of 0.13 pc (first panel on the left of 
Fig. 1). After this, the mass- loss rate and wind velocity are 
changing continuously in the inner boundary giving rise to 
the formation of shocks that develop inside the bow-shock 
cavity formed by the early stationary wind. From then 
onwards the stand-off distance cannot be computed from 
a simple ram pressure balance argument between a free 
streaming stellar wind and the ISM. The ram pressure of 
the ISM is balanced by the ram pressure of the stellar wind 
inside the bow-shock cavity, making it a time dependent 
problem. 

As a consequence of the interaction process a highly 
asymmetric shell is formed that grows in size mainly in 
the downstream direction. The outer bow shock is ra- 
diative for the stellar velocity considered, so the density 
of the swept-up shell is proportional to the square of the 
upstream Mach number and proportional to the ISM den- 
sity. As the mass-loss rate continues (highly modulated 
by the thermal pulses in the star, see Fig. 1 of Villaver et 
al. 2002a) new shells are subsequently formed inside the 
structure generated by the previous stellar wind. During 
most of the AGB evolution (355 000 yr of the - 495 000 
yr that the AGB phase lasts) the leading part of the gas 
is not able to expand far away from the CS. Only the high 
mass-loss rate periods associated with the last two ther- 
mal pulses during the AGB have enough momentum to 
compete efficiently with the previously developed density 
structure. It is at this moment that the leading parts of 
the outer shell growin size. The last two panels show the 
density at the end of the second thermal pulse and in the 
middle of the third one respectively. Because of the rela- 
tive movement, mass is constantly flowing away from the 
head of the bow shock forming cometary structures behind 
the star. 



3.2. The Evolution during the AGB through a Higher 
Density Medium 

In Figure 2 we show the gas evolution for an ISM den- 
sity of n Q = 0.1 cm -3 , lower than the value in the Galactic 
plane (<~ 1 cm -3 ). The snapshots of the evolution have 
been selected at the same times as those shown in Fig- 
ure 1. As expected, the effects of the stellar movement 
are more noticeable, since the ram pressure exerted by the 
ISM depends linearly on the ISM density. The compres- 
sion and deformation of the upstream side of this shell is 
higher than in the case shown in Figure 1. As for the low 
density case, the stagnation radius can be computed from a 
ram pressure balance argument only when the wind is sta- 
tionary in the early evolution (see the first panel of Fig. 2 
where the stagnation radius is a factor 10 -1 / 2 smaller than 
in the first panel of Fig. 1). In the present case, the up- 
stream side of the shell does not increase in size until very 
late in the evolution and the size is always at least 1.5 
times smaller than the one formed when the star evolves 
in a lower density medium. 

3.3. The PN stage 

During the PN phase the stellar wind increases its veloc- 
ity, generating an adiabatic shock that forms the brightest 
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shell of the nebula. In principle, this shell should not be af- 
fected by the interaction with the ISM, unless other mech- 
anisms such as ambipolar diffusion under the presence of 
an external magnetic field are invoked to propagate the 
asymmetry to the innermost regions of the nebula. 

In the right panel of Figure 3 we show the gas density 
when the PN is 2 500 yr old 5 and when the ISM has a den- 
sity of 0.1 cm -3 . The simulated region has been reflected 
around the symmetry axis in order to more clearly show 
the structure. The left panel shows the PN at the same 
age, but when the star is evolving at rest with respect to 
the ISM (see Villaver et al. 2002b). Two bright shells are 
clearly distinguished in both panels of Figure 3. These are 
a spherical innermost shell shaped by the fast wind and 
an external shell which shows a characteristic bow-shock 
configuration in the right panel Figure 3. 

The first effect of the interaction, apart from that on 
the morphology, is that the total size of the outer PN shell 
(halo) is reduced considerably. The radius of the halo in 
the left panel of Figure 3 is 1.8 pc, whereas the outer shell 
of the interacting nebula has a radius of 0.7 pc. The dis- 
placement of the CS with respect to the geometrical center 
of the nebula is 0.1 pc, which corresponds to 20 " at a dis- 
tance of 1 kpc. A less obvious effect is that the formation 
of the outer shell of the PN takes place in a lower density 
environment, since most of the mass ejected during the 
AGB phase is stripped by the ram pressure of the ISM 
and left in the downstream direction of the stellar move- 
ment. The model at rest has 0.7 M of gas above the ISM 
value in the circumstellar envelope (0.43 M that were lost 
by the star and 0.27 M Q of ISM matter swept up by the 
wind), only 0.24 M Q of gas are recovered when the star is 
moving. 

4. DISCUSSION 

We find that the movement of the CS with respect to its 
surrounding medium considerably alters the circumstellar 
gas structure formed during the AGB phase and hence the 
PN formation. Up to now, it has been generally accepted 
that the deceleration of the nebular shell due to its interac- 
tion with the ISM can be more likely observed during the 
late stages of the nebular evolution. We find that when 
the evolution of the stellar wind is properly considered in 
the simulations, the interaction with the ISM cannot be 
studied using simple ram pressure balance arguments and 
plays a major role even during the early AGB evolution. 

From an observational point of view, the presence of 
asymmetries in the haloes of PN that could be related to 
the interaction with the ISM is a common feature 6 . The 
high rate of observed asymmetries cannot be explained 



if all of these PNs arc moving at high velocities through 
high density media. Instead, it can be understood when 
the evolution of the stellar wind is taken into account, pro- 
ducing a much higher efficiency in the interaction. More- 
over, in agreement with our results, Guerrero, Villaver, & 
Manchado (1998) show kinematical evidences that the in- 
teraction of the external shells of PN with the ISM are 
found at all evolutionary stages. 

According to our simulations it may not be necessary 
to invoke the presence of a magnetic field in the ISM to 
explain the asymmetries found in the slow moving (10- 
20 kms" 1 ) PNs Sh 2-216 (Tweedy, Martos, & Noriega- 
Crespo 1995) and NGC 6894 (Zucker & Soker 1997). For 
the model parameters assumed in this paper, no instabili- 
ties appear in the nebula, which is in agreement with the 
analytical results of Dgani & Soker (1994, 1998). 

Another important aspect is that the interaction with 
the ISM considerably reduces the mass of the circumstel- 
lar envelope during the AGB and PN phases due to ram 
pressure stripping. Stellar evolution calculations predict 
that stars with initial masses in the range of <~l-5 M Q 
will end as PN nuclei with masses around 0.6 M Q . Most 
of the mass is lost on the AGB phase and should be easily 
observable as ionized mass during the PN stage. However, 
observations of Galactic PNs reveal on average only 0.2 
M of ionized gas. The mass stripped away by the ISM 
when the star is moving could solve the problem of the 
missing mass in PNs. 

5. SUMMARY 

We have studied the effects the relative movement of the 
CS with respect to the ISM has on the evolution of the 
stellar ejecta during the AGB phase, and on shaping the 
PN. Under the very conservative conditions assumed for 
the ISM density and the CS's velocity (20 kms^ 1 ) we find 
that the PN-ISM interaction provides an adequate mech- 
anism to explain the high rate of observed asymmetries in 
the external shells of PNs. According to our simulations, 
spherical haloes in PNs are expected only if the star is 
at rest with respect to the ISM or if the star is moving 
at low angles with respect to the line of sight, otherwise, 
bow-shock like and cometary structures should be present. 

We are grateful to the referee for very constructive sug- 
gestions. We thank M. L. Norman and the Laboratory for 
Computational Astrophysics for the use of ZEUS-3D. EV 
would like to thank Tariq Shahbaz and Monica Sanchez 
Cubcrcs for fruitful discussions. GGS is partially sup- 
ported by D G APA-UN AM grant IN117799 and CONA- 
CyT. 
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Fig. 1. — The gas density of the structure generated by a star during the AGB phase moving with a velocity of 20 kms~ 1 through an ISM 
with a density of n o =0.01 cm -3 . The intensity scale is logarithmic. All the densities below and above the values shown in the color scale 
have been saturated. The first panel on the left corresponds to 85 000 yr and the subsequent panels are shown at 206 000, 327 000, 387 400 
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Fig. 2. — Same as Figure 1, but for an ISM with a density of 0.1 cm 3 . 
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